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We study ground state properties and dynamics of a dilute ultracold atomic gas in a double well 
potential. The Gaussian barrier separating the two wells derives from the interaction between the 
atoms and a quantized field of a driven Fabry-Perot cavity. Due to intrinsic atom-field nonlinearity, 
several novel phenomena arise being the focus of this work. For the ground state, there is a critical 
pumping amplitude in which the atoms self-organize and the intra cavity field amplitude drastically 
increases. In the dynamical analysis, we show that the Josephson oscillations depend strongly on the 
atomic density and may be greatly suppressed within certain regimes, reminiscent of self-trapping of 
Bose-Einstein condensates in double- well setups. This pseudo self-trapping effect is studied within 
a mean-field treatment valid for large atom numbers. For small numbers of atoms, we consider 
the analogous many-body problem and demonstrate a collapse-revival structure in the Josephson 
oscillations. 

PACS numbers: 42.50.Pq,03.75.Lm,03.75.Nt 
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I. INTRODUCTION 

Outcomes of nonlinearity have extensive consequences 
in various areas of physics leading to phenomena ab- 
sent in their linear counterparts. In quantum mechan- 
ics itplays a crucial role in fields such as nonlinear op- 
tics [l| and ultracold atomic gases [2] . In ultracold atomic 
gases, or Bose-Einstein condensates (BECs), the non- 
linearity stems from interaction among the atoms, and 
cause effects as insulating states of cold atoms in optical 
lattices Q, soliton and vortex formation [4j], collapse- 
revivals of system evolution Q , or self-trappingof atoms 
in optical lattices or double- well (DW) systems \p\ . More- 
over, the inherent phase coherence of the BEG wave- 
function has made it a good candidate for realizing ana- 
logues of the Josephson effect appearing across the bulk 
of two attached superconductors ^TJ . By placing the con- 
densate in a DW potential, the weak tunneling through 
the center barrier brings about oscillations of the BEG 
between the two wells mimicking Josephson oscillation 
in superconductors [8]. Nonlinearity in the BEG DW 
system induces self-trapping and collapse-revivals in the 
Josephson current. Self-trapping emerges for large non- 
linearity and large imbalance of atoms between the two 
wells and manifests itself as blocking of the Josephson 
oscillations. Gollapse-revivals characterize deaths and re- 
births of the oscillations, and become important in the 
quantum regime where fluctuations around the mean- 
field condensate order parameter become relevant. 

In recent years, nonlinearity has turned out to be an 
essential ingredient in two different but highly linked 
matter-light quantum systems, optomechanical cavi- 
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ties and BEG-cavity setups [l^, [HI- Here we focus 
on the second type in which pioneering experiments have 
demonstrated the square root dependence on the number 
of atoms Vljiahi ~ Vn in the vacuum Rabi splitting [lol l. 
Since then, most interest, both theoretically [1^ and ex- 
perimentally [III, [l^ , has indeed been paid to the intrin- 
sic nonlinearity within the system manifested, for exam- 
ple, by optical bistability. In these systems, the cavity 
field acts as an effective potential for the atoms, but the 
atoms, in turn, induce a shift in the index of refraction 
altering the cavity field. As the field adjust accordingly, 
the field cause a back-action on the atoms resulting in 
the nonlinear atom-field interplay. Apart from studies of 
bistability, research has considered also self-organization 
of ultracold atoms in optical resonators, either in setups 
of pumped cavities [13, EB] or pumped atoms [16], as well 
as most recently simulation of the Dicke quantum phase 
transition [iTj . 

In this paper we focus on different outcomes of the 
atom-field nonlinearity compared to earlier works on 
optical bistability. That is, we are not mainly inter- 
ested in multiple-solutions of the equations-of-motion, 
but instead direct our attention on self-organization, self- 
trapping, as well as collapse-revivals in a cavity me- 
diated DW system. The appearance of additional so- 
lutions in nonlinear, in contrast to linear, systems is 
a quite general property. The present work therefore 
aims at more specific outcomes of atom-field nonlinear- 
ity. We notice that BEG DW systems coupled to cav- 
ity fields have been discussed previously [3, [l^, but 
those works, however, did not analyze the situation where 
the cavity field, and thereby the nonlinearity, drives the 
Josephson oscillations. Here self-organization and self- 
trapping are analyzed in the mean-field regime and our 
work goes beyond two-mode approximations commonly 
utilized in DW analysis. We especially show, as for the 
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BEC DW, that the atomic density affects the strength of 
nonhnearity and thereby there is a critical atomic density 
where an effective self-trapping behavior appears. Self- 
organization emerges when the pump amplitude exceeds 
a critical value. At the same instant, the intra cavity field 
amplitude greatly increases similarly to the pump thresh- 
old in the theory of laser. Within the mean- field approach 
we also demonstrate how the system exhibits bistability 
as a result of nonhnearity. Collapse- revival structures de- 
rive from quantum fluctuations, and to tackle such effects 
we consider an effective two-mode many-body model for 
the system which indeed predicts such phenomenon. We 
also point out that the system setup automatically allows 
for quantum nondemolition measurements of the atomic 
dynamics via detection of the output cavity field. 

II. JOSEPHSON OSCILLATIONS AND THE 
MODEL SYSTEM 




FIG. 1: Schematic picture of the traditional DW setup. In 
the symmetric DW, the two onsite energies El and En are 
equal. The left and right solutions are marked by iPl{x) and 
•4}r{x) respectively. 



Before describing the system of the present work, we 
briefly summarize the DW system in order to get a deeper 
understanding for the dynamics analyzed later on in the 
paper. 

In the most simple situation, two wave-functions ijjL{x) 
and SI'S coupled with some strength J. In the sym- 

metric situation, the groundstate and the excited state 
can be written as the symmetric and anti-symmetric so- 
lutions 



(1) 



where Af is the proper normalization coefficient. In 
DW-systems, J represent the tunneling coefficient and 
^l,r{x) are the normalized left well and right well wave- 
functions as depicted in Fig.[T] In general, we can express 
the time-dependent solution vE'(a;,i) = (l)L{t)'4'L{x) + 
4'R{i)i^R{^) leading to a set of two first order coupled 
equations for (jj^ [t) and (j)R{t). The corresponding Hamil- 
tonian is 



HRabi — 



El -J 
-J Er 



(2) 



with El r the onsite energies of the left and right well, 
and J the tunneling strength. Thus, in the symmetric 
well, the detuning 5 — El — Er vanishes. The Rabi 
frequency, characterizing the Josephson oscillations, is 
given by ^Rabi = SV -f 5'^ /Ti, which for the resonant 
well equals 2J /h. Nonzero detunings imply an increased 
Josephson period, as well as causing the amplitude of the 
inversion 
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to decrease from oscillating between -1 and 1 as in the 
case of zero detuning. In deriving we have assumed 
Z{t = 0) = 1, i.e. all atoms initially in the right well. 



A. Cavity mediated double- well for cold atoms 

In both BEC DW experiments of Refs. [l^l, the poten- 
tial barrier separating the left and right wells is obtained 
by dispersive dipole interaction between the atoms and 
an external laser beam. The laser has a transverse Gaus- 
sian mode shape, whose width and amplitude are easily 
adjustable. The large laser intensity makes the light field 
approximately classical. By confining atoms in a res- 
onator, the effective atom-field coupling can be greatly 
enhanced and thereby the dynamics of single atoms can 
be affected by the field even at average photon numbers 
less than unity [2l|. For such low intensities, the light 
field sustained in the resonator must be treated quantum 
mechanically. 

As we are dealing with a coupled bipartite quantum 
system, the state of the atoms will infiuence the state 
of the field and vice versa. Physically, the atomic mat- 
ter wave brings out a change in the index of refraction, 
and as the field adjust to the new index of refraction its 
change will in return affect the atoms leading to an in- 
trinsic nonlinear atom-field interaction. For the present 
system it implies that the effective DW potential directly 
depends on the atomic state. 

The system we have in mind is illustrated in Fig. [2] A 
harmonic potential traps the atoms in all three dimen- 
sions, but the trap frequencies in the y- and 2-directions 
are assumed large enough so that the atoms remain in 
their lowest vibrational states in these directions, and 
we consequently consider motion restricted to the a:- 
direction. Instead of an external laser, we consider the 
field of a cavity to constitute the tunneling barrier. To 
this end we use a Fabry-Perot cavity with its longitudi- 
nal axis along the y-direction. As for a traditional laser 
beam, the Fabry-Perot cavity possesses TiJAfoo-modes 
having transverse Gaussian profiles. The cavity has a 
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FIG. 2: (Color online) System setup. An anisotropic har- 
monic trap confines the atoms such that the motion can be 
considered quasi one-dimensional along the a;-direction. The 
field of a Fabry-Perot cavity, aligned along the y-direction, 
intersects the atomic trap. The TiSMoo-mode shape of the 
cavity induces an effective Gaussian barrier separating the 
harmonic trap into two wells. Furthermore, the cavity is laser- 
driven through one end mirror with an amplitude rj. Cavity 
losses, with decay rate k, are marked by a curly arrow. 



we consider the symmetric DW, and consequently the 
dc (direct current) Josephson effect. The ac (alternating 
current) Josephson effect is assessed by spatially shifting 
either the cavity or the trapping potential. The last non- 
linear term on the right hand side of Eq. (|4]) stems from 
atom-atom collisions and is generally proportional to the 
atomic density and the s-wave scattering length. 

After introducing field losses, k, the corresponding 
Heisenberg-Langevin equations become 



dt 



-(k — iAc)d + rj 



iN / \ii{x)fU{x) dxd + y^dinit), 



2m dx^ 



+ 
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mcu X 
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large enough Q- value to guarantee well separated mode 
frequencies. Nonetheless, cavity losses k are taken into 
account in all our derivations. Only a single cavity mode 
is quasi resonant with the atomic transition under consid- 
eration, and all other modes are therefore ignored. The 
cavity is externally driven with a classical field having 
amplitude rj, and hence, without atoms present the cav- 
ity field would reach a steady state being coherent with 
an amplitude determined by the balancing of losses k and 
pumping rj. The cavity is taken to couple dispersively to 
the atoms and the excited atomic state can hence be 
eliminated adiabatically (22| . In a frame rotating with 
the pump frequency u} „ a nd assuming ultracold atoms, 
the Hamiltonian reads |15| 

Haf = —TiAcd^d — ihrii^d — d^'j 



dxii^x) 



2m dx"^ 



2 2 

mijj X 



(4) 



+u{x)d^d + ^\i!{x)\'' 



*(x). 



where Ac = ujp — uic is the pump-cavity detuning, d) 
(a) the photon creation (annihilation) operator obeying 
standard boson commutation rules, uj the trap frequency, 
and 



U{x) = Woe 



(5) 



is the effective dispersive atom-field coupling with Uq — 
A^/Aa where A is the single photon atom-field coupling, 
Aa = LOp — LOa the pump-atom detuning, and A^r is the 
mode waist. Note that in order to achieve a DW structure 
we restrict the analysis to positive detunings, A^ > 0, but 
interesting effects can also appear for negative detunings 
when the cavity field induces a potential "dimple" on the 
atoms and consequently raises the local phase space den- 
sity. In this work, choosing U{x) centered around x = 



U{x)d''d + gN\^{x)\'' 



Here, din{t) is the Langevin field input noise source 
being (5-correlated; {din{t)din{t')^) = S{t — t'), and 



aLit)a^n{t')) = (a„(t)) 



otherwise 123 
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present work, the infiuence of quantum noise will be as- 
sumed small and hereafter neglected. The effects of such 
fluctuations in a BEG cavity system similar to the one 
we study has been analyzed in Ref . [24J . 

In typical experiments, the characteristic time-scales 
for the field and atoms are substantially different 5(1. 
The field evolution can be assumed to follow the dynam- 
ics of the atoms and it is therefore justified to consider 
the steady state solutions of the field. Explicitly, the 
steady state for the photon number n = a^a is given by 



A, - NY 



(7) 



with the operator Y = J \'^{x)\'^U{x) dx. This equations 
makes clear the modification of the detuning A^ induced 
by the atoms, i.e. the resonance condition Ac = is 
shifted to Ac — NY. Since the atomic field ^{x) is cou- 
pled to the atom number fi, it follows that so is Y and 
then Eq. ([7]) may render multiple solutions of the photon 
number. This same nonlinear effect give rise to quan- 
tum optical bistability fl3, fl3| which can furthermore be 
shown to be analogous to bistable optomechanics in cer- 
tain regimes flTj . 

The intrinsic atom-field nonlinearity is rather different 
from the nonlinearity originating from atom-atom inter- 
actions [2^ . The number of photons enters in the effec- 
tive system potential, and since it depends on the matter 
state the effective potential can be seen as a dynami- 
cal variable. In the BEG DW situation, nonlinearity en- 
ters in the atom-atom interaction term, and within the 
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Thomas- Fermi regime (where kinetic energies can be ne- 
glected) it can be viewed as changing the chemical poten- 
tials in the two wells. Roughly speaking, as an outcome 
in the BEC DW the effective shift due to nonlinearity en- 
ters in the detuning S, while in the cavity DW it modifies 
field strength figs and thereby the tunneling coefficient J. 

Throughout, we try to use realistic experimental pa- 
rameters. More explicitly, we consider the cavity d ecay 
rate k from the experiment of the Esslinger group [l3| . 
and express other rates in terms of this. Lengthscales as 
well as atom numbers are taken as to be experimentally 
realistic. 



III. STATIONARY SOLUTIONS 

In this and the next section we discuss the system at a 
mean-field level, i.e. replacing the atomic operator \E'(x) 
by its mean. The resulting Gross-Pitaevskii equation [26[ 



Q2 _,,2^2 



2m dx^ 



-n,,U{x)+gN\^{x)\'' 



(8) 



VI/ (x). 



for the order parameter can be solved for the minimum 
energy solution by propagating it in imaginary time and 
updating the cavity photon number 



k2 + (Ac - NYf 



(9) 



during the propagation so that it corresponds to the in- 
stantaneous steady-state photon number. Here Y is the 
mean-field counterpart of the operator Y . The corre- 
sponding energy functional is given by [25i [ 



Em 

N 



2to 



dx 



\^{x)\' 



gN 



dx 



kN 



arctan 



N JU{x)\'^/ix)\dx 



In this paper we are not interested in the interaction ef- 
fects due to atomic collisions, but wish to focus on the 
coupling between the cavity mode and the atoms. There- 
fore, we chose g = which can be achieved via Feshbach 
resonances [27| or is approximately valid in a sufficiently 
dilute gas. 

While a direct solution of the Gross-Pitaevskii equa- 
tion is relatively easily found numerically, we found that 
in many instances most of the equilibrium physics (and 
some of the dynamics as well) can also be captured by 



a more transparent variational ansatz. Our ansatz is a 
Gaussian ansatz which can have two peaks, namely 



tp{x) ^ C 



(a:-|-a;o) 



-I- e" 



(11) 



where the prefactor is set by the normalization 
/ dx\il}{x)\'^ — 1 and is given by 



C 



1 



7ri/4 /2a2(l + exp - {^) ) 



(12) 



Substituting this ansatz into the energy functional gives 
us an energy functional E{a, Xq) which can be minimized 
to find the optimal solution. All integrals are sufficiently 
easy to solve analytically, but the resulting expressions 
are too long to be given here explicitly. 

The steady state solutions naturally separate into two 
different regimes, one at the low values of the pump- 
ing strength and the other at large pumping strengths. 
When the pumping strength is small the atomic order pa- 
rameter is nearly a Gaussian centered around the origin. 
In this limit the atomic density overlaps strongly with 
the Gaussian cavity mode function and if the coupling 
strength Uq is quite large, Y can be substantial. Then 
the steady state photon number can be suppressed by 
the TVy-term appearing in the denominator of Eq. Q 
if \NY\ ^ Ac. As the pumping is increased, the pho- 
ton number increases at first roughly cx 77^ followed by a 
pronounced jump in the cavity photon number at some 
critical pump strength. 

This jump coincides with a qualitative change in the 
atomic order parameter as the atomic density splits into 
two separate peaks and develops a minimum at the ori- 
gin. This transition is signaled by a rapid increase in the 
cavity photon number since the overlap of the atomic 
density and the cavity mode function is rapidly reduced 
so that the A^F-contribution in the steady-state photon 
number becomes less important. At the same time the 
effective potential experienced by the atoms becomes a 
stronger double well potential, which pushes the atomic 
density peaks further apart and lowers the overlap with 
the cavity mode function even more. In a way, as one in- 
creases the pumping, one can move from a regime where 
coupling between the atomic order parameter and cavity 
field is strong into a regime where it is weak and the cav- 
ity field mainly acts as an independent optical potential 
acting on the atoms. 

This drastic change in the atomic density of the 
groundstate is an example of self-organization [l^ [28| . 
If the atoms are externally pumped, rather than the cav- 
ity field, the atoms act as a scatter which transfers pho- 
tons from the pump into the resonator. Thereby, there 
is a direct resemblance between pumping of the atoms 
and pumping of the cavity field. The corresponding phe- 
nomenon for pumped atomic systems has been studied 
in great detail and it was found that the transition is de- 
scribed by a second order phase transition ^29j . The same 
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behavior reminds of threshold-pumping of lasers. Pump- 
ing of the laser must exceed a critical value in order for 
lasing to set off. Again, such a critical structure has been 
identified as a second-order phase transition [s^l- More- 
over, it is worth pointing out that in a very recent experi- 
ment, the same kind of transition has been demonstrated 
for a driven lossy cavity in circuit QED (3l| . 

The critical pump strength which separates these two 
self-organized regimes occurs roughly when the barrier 
height becomes larger than the energy scale of the back- 
ground harmonic oscillator i.e. when UggfiUQ ~ Tiui. The 
resonance in the photon number occurs when Ac ~ NY. 
For the Gaussian groundstate of the non-interacting sys- 
tem (which is accurate at small pump strengths) this 
implies NUq/Ac = y^l + cr/A^, where a = ^Jfijmijj is 
the width of the non-interacting (Gaussian) atomic wave- 
function. In this paper we mostly use A^^ = (t/2. 

In Fig. |3] we show an example of the behavior of the 
steady-state photon number as a function of pumping to- 
gether with few examples of the associated atomic order 
parameters as obtained from imaginary time propaga- 
tion of the Gross-Pitaevskii equation. In this figure, we 
also compare the double peaked Gaussian ansatz with the 
results from the Gross-Pitaevskii equation and find the 
agreement to be very good. Using the Gaussian ansatz 
we can, depending on parameters, find two stationary so- 
lutions corresponding to minima of the energy functional. 
The upper branch is the global minimum, while the lower 
one is locally stable. Our Gross-Pitaevskii approach only 
finds the global minimum and in order to find the other 
solution, we would have to impose additional constraints 
on the imaginary time propagation. 

I n ag reement with earlier works on atom-light bistabil- 
ity jl2l.[T3|. the present model supports two or one stable 
values of Ugs- In addition to these solution there exists 
the "middle" branch of the hysteresis curve representing 
an unstable solution corresponding to a local energetic 
maximum. We have not plotted this branch in Fig. [3l 

The above discussion is based on a fairly high value 
[Uq ~ k/100) of atom-field coupling. For a lower value 
of Uq the behavior can be quite different. As an exam- 
ple, we show in Fig. S] results using the double peaked 
Gaussian ansatz and the Gross-Pitaevskii approach with 
otherwise the same parameters, but a much smaller cou- 
pling JJq fa 0.000224k. We chose this value because then 
Ac — NY — and the steady state photon number is on 
resonance for the ideal gas Gaussian groundstate, when 
other parameters are kept the same as in Fig. [31 In this 
case the large overlap between the atomic order parame- 
ter and the cavity mode function acts as to enhance the 
steady state photon number. This is in contrast to the 
behavior at larger couplings where similar large overlap 
implied strong suppression of the cavity photon number. 
In this case we can find no evidence of bistability and the 
double peaked atomic order parameter appears smoothly 
as the pumping strength is increased. At the same time 
the overlap between the cavity mode function and the 
atomic density decreases smoothly. 




FIG. 3: The mean-field photon number as a function of 
r\l K using a double-peaked gaussian ansatz (dashed-line) and 
Gross-Pitaevskii approach (solid line). With some parame- 
ters, the Gaussian ansatz predicts an existence of two station- 
ary solutions corresponding to local minima. Both of these 
solutions are indicated in the figure. As insets we show the 
wavefunction amplitude |^'(a::,i)j using the Gross-Pitaevskii 
equation both at smaller and higher pumping strength. We 
used parameters iV = 10000, «: = 27r x 1.3 MHz, Ac = k, 
l7o = fi:/100, Aa; = O.by'nJrnLj, ui = k/500, and m — 87m. 



IV. MEAN-FIELD DYNAMICS 

In this section we give an overview of the typical fea- 
tures of the system dynamics at the mean-field level. Pre- 
cise details naturally vary based on the used parameters, 
but the basic structures and the physical picture behind 
them is robust with respect to changes in the parameters. 

We can construct wavepackets which are localized ei- 
ther left or right of the barrier by finding the symmet- 
ric (ipsymix)) and anti-symmetric solutions {'ipasymix)) to 
the Gross-Pitaevskii equation. Using these solutions the 
orthogonal solutions localized either to the left or right 
are given by 

IpL^x) = (ipsymix) - asym{x)) / \f2 (13) 



and 



asym 

{x))/V2. (14) 



In Fig. [5] we present an example of the dynamics for an 
initial state which is localized to the right of the barrier. 
Since the barrier height in this case is quite high, the 
overlap between the mode function and the order param- 
eter is small. This remains true even in the course of the 
dynamics and the cavity photon number has only weak 
time-dependence. The atomic order parameter then ex- 
hibits Rabi-fiopping from the right localized state to the 
left localized state with a period which becomes longer as 
the barrier gets higher and/or the cavity mode function 
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FIG. 4: Uss using a double-peaked gaussian ansatz (dashed- 
line) and Gross-Pitaevskii approach (solid line). As insets we 
show the wavefunction amplitude using the Gross- 

Pitaevskii equation both at smaller and higher pumping 
strength. We used parameters A'' — 10000, k — 2n x 1.3 
MHz, Ac = K,Uo = V5k/N « 0.000224 = 0.5^/h/mj, 

u) — k/500, and m — 87u. 



(a) 
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FIG. 5: Dynamics of the atomic order parameter when the 
initial state was localized to the right of the cavity mode bar- 
rier. The upper figure (a) shows the inversion Z{t) together 
with few snapshots of the absolute value |*I'(a:, t)\ of the atomic 
order parameter. The lower plot (b) displays the photon num- 
ber Uss. We used parameters A'' = 10000, k = 2tt x 1.3 MHz, 
77 = 25k, Ac = K, Uo = k/200, A:^ = 0.5^/h/muj, lo = k/500, 
and m = 87u. 



wider. Plotted is the inversion 

Z{t) = 1-2 j \^!{x,t)\^dx, (15) 

and the photon number For the utilized parame- 

ters, t = 100 in dimensionless variables corresponds to 
approximately 6 ms. 

In Fig. [B] we show another example of the dynamics 
for an initial state which is localized to the right of the 
barrier, but with a smaller Uq and larger Ac. In this 
case, the cavity photon number does have pronounced 
time-dependence and the system cannot be described as 
a simple Rabi-flopping anymore. The reason for the dra- 
matically different dynamics is that with the parameters 
used in Fig. [S] the initial barrier height is quite low and 
proportional to 1/(k^ + A^) since the overlap between 
the atomic density and the cavity mode is low. This low 
barrier height enables the wavefunction to populate also 
the barrier region to a greater extent, which increases 
the overlap with the cavity mode. However, such increase 
lowers the (Ac — A^F)^ term in the denominator of Eq. (|9]) 
and causes the photon number to increases toward its 
maximum value oc 1/k^. This increase in turn increases 
the barrier height which tends to drive the atomic order 
parameter away from the region close to the origin. A 
final result of this complicated interplay is a pronounced 
correlated dynamics of the atoms and the cavity field. 
The cavity field reaches a maximum when atoms are, on 
the average, more closely located to the center and is 
a minimum when the order parameter has only a small 
overlap with the cavity mode. 

It should be noticed that while the results in Fig. [S] 



can be analyzed using a two-mode description j6| the re- 
sults presented in Fig. |6] cannot be analyzed in that way 
and a multi-mode description is essential. Since the field 
amplitude is proportional to the photon number, detec- 
tion of the output cavity field intensity would directly 
reveal some properties of the atom dynamics in the cav- 
ity. Moreover, such recording of the output cavity field 
is non-destructive. The idea of utilizing dispersive cav- 
ity interaction for non-demolition measurements has been 
discussed in terms of BEC DW systems [l^ , or for cold 
atoms trapped in optical lattices [32|. These references, 
however, do not consider the cavity field as supplying 
an effective potential and thus is very different from the 
present system. 

When the system was prepared in a localized state 
whose energy is close to the ground state of the double 
well potential, we found a simple Rabi flopping behav- 
ior in Fig. [SJ In Fig. [7] we demonstrate how the state 
can become effectively localized to the one side of the 
double well system when it is prepared in a localized ex- 
cited state. Here we prepared the atoms in a displaced 
Gaussian state and let the system evolve with somewhat 
higher cavity pumping strength than in Fig. [S] As the 
atoms approach the barrier region the photon numbers 
rise, make the barrier higher, and tend to push the atoms 
back. In this case the atoms effectively stay localized 
in one well and the cavity photon number has a regu- 
lar variation which reflects the average position of the 
atoms in a DW system. In this respect, the localiza- 
tion is reminiscent of self-trapping. In regular BEC DW 
systems, increasing the atom number implies enhanced 
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FIG. 6: Dynamics of the atomic order parameter when the 
initial state was a Gaussian locahzed at a;o = 2 to the right 
of the cavity mode barrier and having a width a = 0.8. The 
upper figure (a) shows the inversion Z{t) together with few 
snapshots of the absolute value |^'(a;,t)| of the atomic order 
parameter, while the lower plot gives the photon number Uss- 
We used parameters N = 10000, k = 2-n: x 1. 3 MHz, = 40 k, 
Ac = 3k, Uo = Ac\/5/7V, A^: = 0.5y/h/moj, uj = k/500, and 
m = 87u. 



FIG. 7: Dynamics of the atomic order parameter when the 
initial state was an excited state with a Gaussian profile lo- 
calized to the right of the cavity mode barrier. The upper 
figure (a) shows the inversion Z{t) together with few snap- 
shots of the absolute value of the atomic order pa- 
rameter. Within the time period of this plot, the atoms are 
localized within the right well, i.e. pseudo self-trapped. The 
lower one displays the photon number Uss- We used param- 
eters N = lOOOO, K = 2tt X 1 .3 MHz, r/ = 100k, A^ = k, 
Uo = k/200, A^ = 0.5^h/muj, 



k/500, and m = 87u. 



nonlinearity and a transition between delocalized to lo- 
calized states. This pseudo self-trapping is not perfect; 
the tunneling rate becomes very small but it is still finite, 
while in BEC DW systems it is rather the detuning that 
is affected by the nonlinearity and therefore the oscilla- 
tions are never perfect in this case. Note, however, that 
within typical experimental time scales one would expect 
well established localization. 

In finding the steady state solution to the Gross- 
Pitaevskii equation using the imaginary time propaga- 
tion, we converged to the state corresponding to the up- 
per branch of the Gaussian ansatz and ignored the other 
stationary state. However, even if this state is not the 
lowest energy state it can play an important role in the 
dynamical behavior of the atoms. We illustrate this with 
few examples in Figs. H] and [5] In Fig. [5] we start 
from the steady state at small pumping strengths and 
then ramp it up while in Fig. [5] we do the reverse. Com- 
paring the results indicates a hysteresis type of behavior. 

When the pumping field is ramped up, at first very lit- 
tle happens either for the atomic order parameter or the 
cavity field. There is a broad range of pumping strength 
above the critical pumping strength for the transition to 
a doubled peaked atomic density distribution where the 
system is still dynamically, but not thermodynamically, 
stable. Only when the pumping strength exceeds some 
higher threshold does the system become dynamically 
unstable. When this happens, the atomic order parame- 
ter quickly splits, the cavity photon number is suddenly 
increased, and the energy of the atoms is increased con- 



siderably in this abrupt process. 

This behavior is easy to understand, since here we start 
with a single peaked atomic order parameter which has a 
substantial overlap with the cavity mode function. This 
means that the photon numbers, and hence the barrier 
height, are strongly suppressed by the iVF-term in the 
denominator. Only once the pumping strength becomes 
so strong that the barrier height is increased sufficiently 
to split the wavefunction, does the picture change. Then 
the overlap parameter Y is quickly reduced and the pho- 
ton number (and the barrier height) increase. This in 
turn splits the atomic order parameter even more and 
results in the abrupt dynamics we observe. 

On the other hand, when the field is ramped down (Fig. 
ini) the initial order parameter has only small overlap 
with the cavity mode function. Reducing the pumping 
strength will then reduce the barrier height smoothly al- 
lowing the order parameter to contracts toward the single 
peaked order parameter more smoothly. Eventually the 
overlap between the order parameter and the cavity mode 
becomes so strong that the NY-teim in the denomina- 
tor of Eq. ^ reduces the photon number more strongly 
than the simple reduction of the pumping strength would 
suggest. 

Up till now we considered an atom number N — 100 
which is of the correct order of magnitude with present 
experiments [To| . The thermodynamical limit is given 
by letting N and the effective mode volume V tend to 
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(a) (a) 




t[1/(nl t[1/to] 



FIG. 8: The pump strength (a), mean-field photon num- 
ber (b), and snapshots (insets) of the atomic order param- 
eter |»I'(a:;,t)| as a function of time when the pump strength 
is sweeped from low to high values. After a certain pump 
strength, the barrier separating the left and right wells be- 
comes sufficiently large to split the atom density and simul- 
taneously the photon number makes a drastic change. We 
used parameters TV = 10000, k = 2-k x 1.3 MHz, Ac = k, 
Uo = k/200, Aa; = 0.5y^h/muj, uj = k/500, and m = 87u. 



FIG. 9; The pump strength (a), mean-field photon number 
(b), and snapshots (insets) of the atomic order parameter 
\'i'{x,t)\ as a function of time when the pump strength is 
sweeped from high to low values. In comparison to Fig. [8l 
the transition is here much smoother. We used the same 
parameters as in Fig. [8] i.e. N = 1 0000, k = 2n x 1.3 MHz, 
Ac = K, Uo = k/200. Ax = 0.5^/h/rruu, uj = k/500, and 
m — 87u. 



infinity while keeping the density p — N/V fixed. In 
the driven cavity QED system one has Uo oc and 
T] oc VV [s^l, and it follows that the effective potential 
nssU{x) depends solely on p, not N and V independently. 
Consequently, the same characteristics is obtainable for 
other atom numbers N by suitably choosing eta and Uq. 



DYNAMICS BEYOND MEAN-FIELD 



Since it is known that often the relevant physics can 
be captured by a two-mode model we consider such an 
approach in this section which in particular allows us 
to explore physics beyond mean-field. Especially for 
deep enough DWs, it is legitimate to assume that the 
ground state can be written as in Eq. ([T]), i.e. 5*0 (a^) = 
^ {iPl{x) + iPr{x)) with i^LMi^) localized states in the 
left and right wells. As in the previous sections, the left 
and right functions 'ipL,Rix) are taken to be Gaussian in 
form. We pick the center of these to coincide with the 
minima of the effective potential 



The width, on the other hand, is obtained from mini- 
mizing the corresponding energy functional. In princi- 
ple, xo could be kept a variational parameter, but in or- 
der to simplify the estimations of the parameters of our 
two-mode model, we impose this additional assumption. 
Relaxing this constraint, is straightforward, but is not 
expected to change our results in an important way, es- 
pecially when the barrier height is large compared to the 
energy scale of the harmonic trap. 



givmg 



Xo 



A:,Whi 



U{x)ns 



I'm/or 



(16) 



(17) 



Thus, the width a is found by minimizing the sys- 
tem's energy functional. This results in a set of coupled 
nonlinear equations for a and Uss which must be solved 
self-consistently. In particular, multiple solutions of the 
equations may exist in certain parameter regimes. Intro- 
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ducing the overlap integrals 

-2 



Eo 
El 



(18) 



^0 



Si 



\il)i{x)\ X dx 



■>Pl{x)'ipR{x)x'^dx, 



the groundstate energy within the present ansatz takes 
the form I25| 



E 
N 



= En + Ei+Sn + Si 



r , ('A,-UoN{Jo + Ji) 

circtcLii 

kN 



K 



(19) 



We artificially impose orthogonality of the left and right 
well wave-functions (J ip'^(x)ipR{x)dx — 0) in order to 
avoid unphysical contributions, and then get 



En = 



Jn = 



A., 



VAT 



(T 

~2 



El = 
J I = 

Si - 



A., 



n/aT 



(20) 

As pointed out, in order to analyze effects beyond mean- 
field we assume the system parameters to be such that 
we can impose a two-mode approximation. Thus, we ex- 
pand the atomic operators as ^'(a;) = bLipL{x)+bRipRix), 
where (^1, ji) annihilates (creates) an atom in well 
L, R, and 4'l,r{x) are determined as above. The two 
functions '>pL,R.i^) are operator valued since they depend 
on the cavity field amplitude. As an outcome, in deriv- 
ing equations-of-motion or an effective Hamiltonian for 
the atoms when the cavity field is eliminated, one has 
to take ordering between non-commuting operators into 
account [l^. To do so, we will assume Jq ^ Ji and leave 
out the Ji cross-term in Y for the steady-state photon 
number of Eq. ([T]). For the present potential U{x), this 
is not always justified but nevertheless it holds in large 
parameter regimes. In this work we will not present the 



full derivation of the effective Hamiltonian, but refer to 
Ref. [15^ for details. To order we find a second 

quantized Hamiltonian for the atoms 



Hbh = {Eo + So)N + fiN) - t{N)B, 



(21) 



where 



N = fiR + fiL ^ btbR + blbL, 



B = blbL + blbn, 
t{N) = -El - Si - 



fiN) = — arctan 

K 



WUoJi 



+ (A, - UoJqNY 
( A,-UoJoN\ 



(22) 



and the coefficients are given in Eq. (1201 
Transforming the operators as 



'h' 


1 


■ 1 1 ■ 




'k' 




"71 


1 -1 







(23) 

"J L J 
the Hamiltonian is diagonalized 

{Eo + So)N + f{N)-t{N){n+-h^), (24) 



H 



BH 



with the new number operators n+ = ^+^-(- and ri_ = 
6lS_. For given atom number N ^ there are A^-|-l equidis- 
tant energy levels separated by t{N). In this case we 
recover perfect Josephson oscillations with an oscillation 
frequency Vl.Rabi{N) = 2\t{N)\/fi. For atomic states with 
an uncertain number of atoms, the various A'^'s will in- 
duce a collapse in the Josephson oscillations as the con- 
tributing terms move out of phase. Such a collapse was 
discussed for the regular EEC DW in Ref. @, where it 
derives from the atom-atom interaction term. For mod- 
erate or large atom numbers it is appropriate to assume 
an initial coherent atomic state 



IV') 



, . Ajri/2 



(25) 



Here, A" = (A^) is the average number of atoms and the 
state jn,™) gives the number n of atoms in the left well 
and m atoms in the right well. For N ^ oo, the relative 
uncertainty SN = AN/N, where AA^ = {{N - N)^), 
goes to zero representing the "classical" (mean- field) 
limit. Depending on the particular A^-dependence of 
t{N), the different Josephson oscillation terms may re- 
turn in phase causing the system to revive The 
revival time Tr can be estimated as the time it takes 
for consecutive terms to build up a 2tt phase difference 
Tr[ioj{N + 1)- ujj{N)] = 27r giving 



N=N 



(26) 
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The A'^-dependence of t{N) is supposedly weak in the 
large atom limit implying — >■ oo as iV — >■ oo, as ex- 
pected in the classical limit. Collapse-revival patterns are 
a widespread phenomena in physics and have especially 
been studied in the vibrational dynamics of molecules |34l | 
and in cavity QED j33l |. as direct proofs of quantization 
of either the molecular vibrations or of the electromag- 
netic field. In the BEC DW system the collapse-revivals 
derive from the squared atom number operators, n| ^, 
while in standard cavity QED it is typically an outcome 
of a square-root dependence of photon numbers in the 
Jaynes-Cummings model (ssj , Vh. Here, the atom num- 
ber dependence of t{N) is presumably more complex and 
one thereby expects a less pronounced collapse-revival 
structure. 



(a) 




(b) 



t [l/a] 
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FIG. 10: Time evolution of the inversion Z{t) for an initial 
coherent state in the right well with N = 50 (a) and iV — 100 
(b). The common parameters are, k = 27r x 1.3 MHz, Ac — k, 
•q = 2k, Uo = k/5. Ax ~ 0.1 ^m, and uj = k/500. 



In Fig. 1101 we display the time evolution of the many- 
body inversion 



-■MB 



{ur) - (nii) 
N 



(27) 



for an initial coherent atomic state in the right well with 
TV = 50 (a) and N = 100 (b). At short times, t < 0.3, 
a few Josephson oscillations persists before the collapse. 
The collapse time is approximately the same for both 
examples, which is a general property [33| . The first 
collapse period lasts until i ^ 2.7 for iV = 50 and until 
t ~ 3.2 for N = 100. The iV-dependence in the revival 
times is as well a general feature as argued above; large 
N values imply long revival times [33|. In the upper 
plot (a), a sequence of collapse- revivals appear, however 
less clear as time progresses. For the lower plot (b), on 
the other hand, only a single collapse period is visible. 
When the time-span of revivals begin to overlap, which 
start to happen after roughly 1 ms in (a) and already 
after the first revival in (b), super-revivals may occur 



due to higher order interferences [3g]. Slight signatures 
of super-revivals can be seen in (b) in the modulated 
oscillations after t ^ 10. We note that for the present 
results the time-spans are of the same order as those for 
single Josephson oscillations in Fig. [TOl this derives from 
the much weaker pump amplitude r] in these examples 
giving a lower barrier between the two wells, i.e. shorter 
tunneling times. 

As pointed out in the mean-field section, we have ne- 
glected effects arising from atom-atom nonlinearities by 
letting g = 0. The quadratic terms in atom numbers are 
known to render collapse-revivals f34| , and with both 
nonlinearities present (atom-atom and atom-field inter- 
actions) one would see a competition between the two 
mechanisms. In general, revivals become less frequent 
since they are only possible when both effects simultane- 
ously support revivals. We have verified these conclusions 
numerically by including atom- atom interactions into our 
calculations. 

The analysis of this section relies on a static assump- 
tion, i.e. the effective potential is considered indepen- 
dent on time. The collapse-revival structure is solely 
an outcome of the uncertainty of atom number and not 
of the instantaneous state of the atoms. The full time- 
dependent many-body problem is certainly interesting, 
but in this work we focus on general novel phenomena 
inherent in the atom-field nonlinearity. Even though the 
results of these section have been derived within some 
assumptions, we believe that the general structure sur- 
vives also in more rigorous analyzes. Such approaches 
would for example reveal how the steady-state photon 
number n^s — (fiss) evolves in time. Nevertheless, it is 
believed that Uss encodes the properties of the atomic 
evolutions exactly as it did in the previous section study- 
ing mean-field dynamics. That is, one would expectedly 
find a similar collapse-revival pattern in Uss as the one 
found in the atomic inversion. 

It is important to note that the collapse-revival phe- 
nomenon is a pure quantum effect. By adding fluctua- 
tions around the means ^'(x) of the previous two sections 
one would typically encounter a collapse in the oscilla- 
tions. However, the revivals depend strongly on the par- 
ticular A^-dependence of t{N) and therefore random fluc- 
tuations around 4'(x) would not predict revivals. In other 
words, capturing the revival structure from some per- 
turbed mean-field approach is only possible if one could 
carefully chose the perturbing fluctuations. 



VI. CONCLUSIONS 

Nonlinear Josephson oscillations have been studied in 
a DW system of ultracold bosonic atoms. Contrary to 
regular nonlinearity, arising from atom-atom interaction 
in these types of systems, the nonlinearity we considered 
derives from intrinsic interaction between the atoms and 
a quantized cavity field. In particular, we demonstrated 
the appearance of pseudo self-trapping, self-organization. 
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as well as collapse-revivals. In other words, the present 
work shows that these phenomena are rather general and 
not restricted to only certain kinds of nonlinearities. 

The system parameters were chosen in agreement 
with current experiments indicating that the phenom- 
ena should be experimentally realizable. We furthermore 
showed or argued that the output cavity field, propor- 
tional to Uss, provides a direct handle of the atomic dy- 
namics. Since the output field is regarded as losses, mea- 
surement of it is nondemolition by nature. Hence, the 
various evolution regimes can be traced down without 
standard destructive measurements such as time-of-flight 



or fluorescence detection. 
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